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Abstract: This paper is devoted to the study of the asymptotic dynamics of 
a class of coupled second order oscillators driven by white noises. It is shown 
that any system of such coupled oscillators with positive damping and coupling 
coefficients possesses a global random attractor. Moreover, when the damping 
and the coupling coefficients are sufficiently large, the global random attractor 
is a one-dimensional random horizontal curve regardless of the strength of 
the noises, and the system has a rotation number, which implies that the 
oscillators in the system tend to oscillate with the same frequency eventually 
and therefore the so called frequency locking is successful. The results obtained 
in this paper generalize many existing results on the asymptotic dynamics for 
a single second order noisy oscillator to systems of coupled second order noisy 
oscillators. They show that coupled damped second order oscillators with 
large damping have similar asymptotic dynamics as the limiting coupled first 
order oscillators as the damping goes to infinite and also that coupled damped 
second order oscillators have similar asymptotic dynamics as their proper space 
continuous counterparts, which are of great practical importance. 
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1 Introduction 

This paper is devoted to the study of the asymptotic dynamics of the following system of second 
order oscillators driven by additive noises: 



where j E Z% := {j = . . . , j d ) G Z d : 1 < ji, ■ ■ ■ , jd < N}, Uj is a scalar unknown function of 
t and u = (ui, it2, • • • , u N d) T , a and K are positive constants, A is an N d x N d matrix and (Au)j 
stands for the jth component of the vector Au, /3 6 R, g is a periodic function, fj and €j are 
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constants, and {Wj(t)}j &Z d are independent two-sided real- valued Wiener processes. Moreover, 
A and g satisfy 

(HA) A is an N d x N d nonnegative definite symmetric matrix with eigenvalues denoted by X{, 
i = 0, 1, . . . , N d — 1 satisfying that 

= Ao < Ai < • • • < X N d_ 1 , 

Ao is algebraically simple, and (1, . . . , 1) T £ M. Nd is an eigenvector corresponding to Xq = 0. 
(HG) g G C 1 (M,1R) has the following properties 

g(x + K)=g(x), \g{x)\<ci, \g'{x)\ < c 2 , Vx € R, 

where c\ > 0, C2 > and k > is f/ie smallest positive period of g. 

System (II. ip appears in many applied problems including Josephson junction arrays and 
coupled pendula (see [11], [13], [22], [27], etc.). Physically, a in (jl.ip represents the damping of 
the system and K is the coupling coefficient of the system. (11.11) then represents a system of 
N d coupled damped oscillators independently driven by white noises. 

System (jl.ip also arises from various spatial discretizations of certain damped hyperbolic 
partial differential equations. For example, the -/V d x N d matrix A in (jl.ip includes the dis- 
cretization of negative Laplace operator —A with Neumann or periodic boundary conditions 
defined as follows: 

(Au)j = (Au)(j 1)i2) ... )jd) = j^[2duj - ^(jx+i.^vj'd) ~ u (juh+h-,jd) u Ui,h,-,j d +^ 

~ u {h-hh,-,id) ~ u Ui,h-l,-,jd) — ■ ■ ■ — u 0'i.J2,-.,id-l)] ' 
with Neumann boundary condition 

U (jli— Ji-l.OJi+Xi— iJd) = U (jli— iJi-lAiii+X)— J'd)' 

u (ji,...,ji-i,N+l,j i+ i,...,jd) = u (h,-di-i,N,ji+i,-,jd) 
or periodic boundary condition 

u {ji,—,3i-i,0,ji+i,-,3d) ~ U {h,-,3i-uN,ji + i,...,j d )i 

U (jl,—ji-l,N+l,ji + i,...,j d ) ~ U {h,—Ji-lA>3i+l,-Jd) 

for j = (ji,...,j d ) e and i = l,...,d. Thus, ([□]) with Uj = u{jih,--- ,jih,--- ,j d h) 
{h = L/N), A being as above, fj = /, €j = e, and Wj = W is a spatial discretization of the 
following problem 

dii + adu - KAudt + (3g{u)dt = fdt + edW, in U x R + (1.2) 
with Neumann boundary condition or periodic boundary condition, i.e., 

— = ondU xM + 
on 



or 

du 



M r = U r., 



du 

dx-; 



dxj 



, j = l,...,d, 

F j + d 
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where Tj = dU n {xj = 0}, = dU n {scj = L}, j = 1,. . . ,d and 17 = nf=i(°> L )- Note that 
if g(u) = sinn, (11.21) is the so called damped sine-Gordon equation, which is used to model, for 
instance, the dynamics of a continuous family of junctions (see |25j). 

Two of the main dynamical aspects about coupled oscillators and damped wave equations 
considered in the literature are the existence and structure of global attractors and the phe- 
nomenon of frequency locking. A large amount of research has been carried out toward these 
two aspects for a variety of systems related to (jl.ip . See for example, [El [18], [20], [22] [28] for 
the study of coupled oscillators with constant or periodic external forces; |12j [T9] [2T] [25] [26] 
for the study of the deterministic damped sine-Gordon equation; [8] HE] [23] for the study of 
coupled oscillators driven by white noises; and [101 [Ml [29] for the study of stochastic damped 
sine-Gordon equation. Many of the existing works focus on the existence of global attractors and 
the estimate of the dimension of the global attractors. In [8, 23] [23], the existence and structure 
of random attractors of stochastic oscillators and stochastic damped wave equations are studied. 
In particular, the asymptotic dynamics of a single second order noisy oscillator, i.e., (jl.ip with 
N = 1, is studied in [23]. The author of [23] proved the existence of a random attractor which is 
a family of horizontal curves and the existence of a rotation number which implies the frequency 
locking. In [8], the authors considered a class of coupled first order oscillators driven by white 
noises. Among those, the existence of a one-dimensional random attractor and the existence of 
a rotation number are proved in [8] . The system of coupled first order oscillators considered in 
[8] is of the form 

duj + K(Au)jdt + (3g(uj)dt = fjdt + edWj, j £ Z%. (1.3) 
Note that, by resealing the time variable by t — > — , (jl.ip becomes 

-duj + duj + K(Au)jdt + f3g(uj)dt = fjdt + edWj, j G Z%. (1.4) 

Hence, (jl.3p can be formally viewed as the limiting system of (jl.ip as the damping coefficient a 
goes to infinite. In [23]; the authors investigated the existence and structure of random attractors 
of damped sine-Gordon equations of the form (jl.2p with Neumann boundary condition, which 
is a space continuous counterpart of (jl.ip as mentioned above. 

However, many important dynamical aspects including the existence of global attractor and 
the occurrence of frequency locking have been hardly studied for coupled second order oscillators 
of the form (jl.ip driven by white noises. It is of great interest to investigate the extent to 
which the existing results on asymptotic dynamics of a single second order noisy oscillator may 
be generalized to systems of coupled second order noisy oscillators. Thanks to the relations 
between (jl.ip and (|1.3|) and between (jl.ip and (|1.2p . it is also of great interest to explore the 
similarity and difference between the dynamics of coupled damped second order oscillators and 
its limiting coupled first order oscillators as the damping coefficient goes to infinite and between 
the dynamics of coupled damped second order oscillators and their proper space continuous 
counterparts. The objective of this paper is to carry out a study along this line. In particular, 
we study the asymptotic or global dynamics of (jl.ip . including the existence and structure of 
global attractor in proper phase space and the success of frequency locking. 

In order to do so, as usual, we first change (jl.ip to some system of coupled first order random 
equations. Assume N >2 and d > 1 (N = 1 reduces to the single noisy oscillator case considered 
in [23]). Let u = («i) i6Z * , g(u) = (^j))^, / = {Ji)j&.%> = (^^(^iez^- Then > 
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(jl.ip can be written as the following matrix form, 



du + adu + KAudt + (3g(u)dt = fdt + dW(t). (1.5) 

Let 

Uj = n = {uj e C(R,R) : w (0) = 0} 

equipped with the compact open topology, Tj = £>(f^o) be the Borel cr-algebra of f^o and Pj be 
the corresponding Wiener measure for j £ Z^. Let = Ylj& d %' De * ne product <7-algebra 
on f2 and P be the induced product Wiener measure. Define (6t)t<=m. on Q via 

e t uj{-) =u{- + t) - uj(t), tel. 

Then, (f2, F, P, (St)^^) is an ergodic metric dynamical system (see QJ). Consider the Ornstein- 
Uhlenbeck equation, 

dz + zdt = dW(t), z£R N \ (1.6) 

Let z(6tuj) = (zj(9tUj))j eZ d be the unique stationary solution of (|1.6|) (see [TJ [5J [U] for the 
existence and various properties of z{-)). Let v = u — z{0toS). We obtain the following equivalent 
system of (jl.5p . 

| v = -KAu -atv + f- j3g{u) + (1 - a)z(6»^). 

To study the global dynamics of (jl.ip , it is therefore equivalent to study the global dynamics 
of (|1.7p . Observe that the natural phase space for flT7]) is E := 1^ xt" with the standard 
Euclidean norm. Thanks to the presence of the damping, it is expected that (ll.7p possesses a 
global attractor in certain sense. However, due to the uncontrolled component of the solutions 
along the direction of the eigenvectors of the linear operator in the right of (jl.7|) corresponding 
to the zero eigenvalue, there is no bounded attracting sets in E with the standard Euclidean 
norm, which will lead to nontrivial dynamics. There is also some additional difficulty if one 
studies (jl.7p in E with the standard Euclidean norm due to the zero limit of some eigenvalues 
of the linear operator in the right of (11. 7p as a — > oo. The later difficulty does not appear for 
coupled first order oscillators studied in [8] and for a single noisy oscillator considered in [23] , 
We will overcome the difficulty by using some equivalent norm on E and considering (II. 7p in 
some proper quotient space of E and prove the existence of a global random attractor as well 
as the existence of a rotation number of (|1.7p . 

To be more precise, let 

C ={-KA -al)- < L8 > 
By simple matrix analysis, the eigenvalues of C are given by (see [14} [20] for example) 

»f= - a±Va : < = 0,1 JV--!. (1.9) 



Note that fi^ = 0, which requires some special consideration for the solutions along the direction 
of the eigenvector t]q = (1, . . . , 1, 0, . . . , 0) T corresponding to /i^j". We overcome this difficulty by 
considering (jl.7p in the cylindrical space E\jnr]^L x E2, where E\ = spanj^o}, E2 is the space 
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spanned by all the eigenvectors corresponding to non-zero eigenvalues of C (see section 4 for 
details). We then prove 

(1) For any a > and K > 0, system (jl.ip possesses a global random attractor (which is 
unbounded along the one- dimensional space E\ and bounded along the one-codimensional space 
E2) (see Theorem \4-£\ Corollary \4-3\ and Remark 4-4)- 

It is expected physically that when the damping coefficient a — > 00, the dynamics of (11 . 71) 
becomes simpler or the structure of the global attractor of (|1.7p becomes simpler. However, 
pf — > as a — > 00 for i = 1, 2, • • • , N d — 1, which gives rise to some difficulty for studying 
the structure of the global attractor in E with the standard Euclidean norm. We introduce 
an equivalent norm on E to overcome this difficulty (see section 3 for the introduction of the 
equivalent norm, the choice of such equivalent norm was first discovered in [15]) and prove 

(2) When a and K are sufficiently large, the global random attractor of (jl.ip is a one- dimensional 
random horizontal curve (see Theorem \5.3\ and Corollary \5.4\ ), and the rotation number (see Def- 
inition [6A\) of (jl.ip exists (see Theorem 1 6. 31 and Corollary \6.4\ )- 

Note that roughly a real number p £ M is called the rotation number of (jl.ip or (11. 7j) if for 
any solution {uj(t)}j e %d of (jl.ip . the limit lim^oo exists almost surely for any j £ and 

lim ^2^1 = p for a.e. uj G VL and j = 1, 2, • • • , N d 

>oo t 

(see Definition 16.11 and the remark after Definition 16. ip . Hence if (jl.ip has a rotation number, 
then the oscillators in the system tend to oscillate with the same frequency eventually and 
therefore the so called frequency locking is successful. 

(1) and (2) above are the main results of the paper. They make an important contribution 
to the understanding of coupled second order oscillators driven by noises. Property (1) shows 
that system (jl.ip is dissipative along the one-codimensional space E2. By property (2), the 
asymptotic dynamics of (jl.ip with sufficiently large a and K is one dimensional regardless of 
the strength of noise. Property (2) also shows that all the solutions of (jl.ip tend to oscillate 
with the same frequency eventually almost surely and hence frequency locking is successful in 
(jl.ip provided that a and K are sufficiently large. 

The results obtained in this paper generalize many existing results on the asymptotic dynam- 
ics for a single damped noisy oscillator to systems of coupled damped noisy oscillators. They 
show that coupled damped second order oscillators with large damping have similar asymptotic 
dynamics as the limiting coupled first order oscillators as the damping goes to infinite and hence 
one may use coupled first order oscillators to analyze qualitative properties of coupled second 
order oscillators with large damping, which is of great practical importance. They also show 
that coupled damped second order oscillators have similar asymptotic dynamics as their proper 
space continuous counterparts and hence one may use finitely many coupled oscillators to study 
qualitative properties of damped wave equations, which is of great practical importance too. 

The rest of the paper is organized as follows. In section 2, we present some basic concepts and 
properties for general random dynamical systems. In section 3, we provide some basic settings 
about (jl.ip and show that it generates a random dynamical system. We prove in section 4 the 
existence of a global random attractor of the random dynamical system <ft generated by (jl.ip for 
any a > and K > 0. We show in section 5 that the global random attractor of <j> is a random 
horizontal curve and show in section 6 that (jl.ip has a rotation number, respectively, provided 
that a and K are sufficiently large. 
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2 Random Dynamical Systems 

In this section, we collect some basic knowledge about general random dynamical system (see 
[TJ 0] for details). Let (X,d) be a complete and separable metric space with Borel u-algebra 
B(X). 

Definition 2.1. A continuous random dynamical system over (f2, J 7 , P, (0t)teM.) is a (B(R + ) x 
T x B(X),B(X)) -measurable mapping 

<p : M. + x S7 x X — > X, (t, lj, x) \-¥ (p(t,u,x) 

such that the following properties hold: 

(1) ip(0,u),x) = x for all uj G Q; 

(2) tp(t + s, uj, •) = (p(t, 9 s uj, tp(s, uj, •)) for all s, t > and uj G S7; 
f5j (p(t,uj,x) is continuous in x for every t > and we!!. 

For given x £ X and £7, F C X, we define 

d(x, F) = inf d(x, y) 

y&F 

and 

d H {E,F) = sup d(x,F). 

x&E 

dfj(E, F) is called the Hausdorff semi-distance from E to -F. 

Definition 2.2. (1) A set-valued mapping uj h- > Z?(u;) : — >• 2^ zs said to 6e a random set 
if the mapping ui \- )■ d(x, D(uj)) is measurable for any x £ X. If uj i— >■ d{x,D{ui)) is 
measurable for any x £ X and D(uf) is closed (compact) for each uj £ ft, then uj h-» -D(w) 
is called a random closed (compact) set. A random set uj \-t D(u>) is said to be bounded 
if there exist xq £ X and a random variable R(uj) > such that 

D(uj) C {x £ X : d(x,x ) < R(uj)} for all uj G SX 

(2) A random set ui \— > D(u) is called tempered provided that for some xo £ X and P-o.e. 
a; G SI, 

lim e"' 3 * sup{d(&, x ) : 6 G £>(0_tw)} = for all > 0. 

i— >oo 

f5j ^4 random set uj \— > B{ui) is said to be a random absorbing set if for any tempered random 
set uj i — y D(oj), there exists to(oj) such that 

tp(t,0- t u,D(9- t u)) C B(oj) for all t > t (uj), ui £ Q. 

(4) A random set oj \— > R\{oj) is said to be a random attracting set if for any tempered random 
set uj i — y D(uj), we have 

lim dH(<p(t, D(9-tUj), B\(oj)) = for all uj £ O. 

t— 'too 

(5) A random compact set uj \- > A(uj) is said to be a global random attractor if it is a random 
attracting set and (p(t,uj, A{ui)) = A(8tuj) for all uj £ £1 and t > 0. 
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Theorem 2.3. Let ip be a continuous random dynamical system over (Q, J 7 , P, (9t)teu)- If there 
is a random compact attracting set oj \— > B{oS) ofcp, then lo \— > A(uS) is a global random attractor 
of ip, where 

t>0r>t 

Proof. See [Hi]. □ 

3 Basic Settings 

In this section, we give some basic settings about (jl.ip and show that it generates a random 
dynamical system. 

First, let Y = {u,v) T and F(O t u,Y) = (z(O t u),f- Pg{u) + (1 - a)z{9 t uj)) T . System (fL7]) 
can then be written as 

Y = CY + F(9 t oo,Y), (3.1) 

where C is as in (|1.8p . 

Recall that z{6toS) = (zj(6 t uj)) je ^d is the unique stationary solution of (11.60 . Note that the 
random variable is tempered and the mapping t i-> Zj{6 t oo) is P-a.s. continuous (see 

[H [2]). More precisely, there is a ^-invariant f2 C fi with P(f2) = 1 such that t \— > Zj(0fUj) is 
continuous for u; G f2 and j G Z^-. We will consider (jl.7p or (|3.ip for w £ fi and write O as $7 
from now on. 

Let E = R Nd x R Nd and F w (r,Y~) := F(0 t w,Y'), then F w (-,-) : R x E —?■ E is continuous 
in t and globally Lipschitz continuous in Y for each uj G O,. By classical theory of ordinary 
differential equations concerning existence and uniqueness of solutions, for each oj G and any 
Yq G E, (13. ip has a uniqueness solution Y(t, u),Yq), t>0, satisfying 

Y(t,uj,Y ) = e ct Y + [ e cit ' s) F(e s uj,Y(s,uj,Yo))ds, t > 0. (3.2) 
Jo 

Moreover, it follows from [1] that Y(t, uj, Yq) is measurable in (t, uj, Yq). Hence (|3,ip generates a 
continuous random dynamical system on E, 

Y:R + xQxE^E, (t,u,Y ) t->Y(t,u,Y ). (3.3) 

Define a mapping : R + x $7 x E — > E by 

<j>[t,u,<M = Y(t,oj,Y (uj)) + (0,z(6 t uj)) T , (3.4) 

where 4>q = (uo,ui) T G E 1 and lo( w ) = («0)^i — z i UJ )) T ■ Then is a continuous random 
dynamical system associated with the problem (jl.lj) on 

Recall that the eigenvalues of C are given by (see |144 [20] for example) 

,t= ~ a± ^- AK \ i = 0,l,...,N«-l. (3.5) 

By (|3.5p . C has at least two real eigenvalues and —a with eigenvalues 770 = (1, ...,1,0,..., 0) T , 
77-1 = (1, . . . , 1, — q, . . . , — a) T G E, respectively. Let E\ = span{r/o}, -E-i = span{?]_i}, i?n = 
E\ + and E22 = ^ij, the orthogonal complement space of E\\ in E, then E 1 = £n © #22 ■ 
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To control the unboundedness of solutions in the direction of rjo, we will study (|3.ip in the 
cylindrical space Ei/kt]qZ x E 2 , where E2 = E-\ © £"22 (see Section 4 for details). 

Observe that the Lipschitz constant of F with respect to Y in E with the standard Euclidean 
norm is independent of a > 0. But [if — > as a — > 00 for i > 1, which gives rise to some 
difficulty for the investigation of (|3.ip in E with the standard Euclidean norm. To overcome the 
difficulty, we introduce a new norm which is equivalent to the standard Euclidean norm on E. 
Here, we only collect some results about the new norm (see [15} 120] for details). 

Define two bilinear forms on En and E22, respectively. For Yi = (ui,Vi) T G -En, i = 1,2, let 

2 

or ol ol 

(Yi, Y 2 )e u = -^(ui,u 2 ) + (-ui +vi, -u 2 +v 2 ), (3.6) 

where (•, •) denotes the inner product on R 7 ^, and for Y{ = (ui,Vi) T G £22, i = 1,2, let 

2 

or ol ol 

(^1,^2)^32 = (KAu 1 ,u 2 ) + (— - 5KX 1 )(u 1 ,u 2 ) + (-1*1 + vi, -u 2 + v 2 ), (3.7) 

where 5 G (0, 1]. It is easy to check that the Poincare-type inequality 

(Au, u) > Ai ||-u|| 2 , V Y = (u, v) T G E22 

holds (see |2U] for example), where || • || is the standard Euclidean norm on R^. Thus (|3.7p is 
positive d< 
we define 



positive definite. For any Yi = +y} 2) G E, i = 1, 2, where y} 1} , Y 2 (1) G E n , y} 2) , Y 2 (2) G E 22 



(Yx,Y 2 ) E = (Y^, Y 2 (1) >£u + (y} 2 \yP) E22 . (3.8) 
Lemma 3.1 ( [20] ) . (1) (I3.6P and (|3.7p define inner products on En and E22, respectively. 

(2) (13. 8p defines an inner product on E, and the corresponding norm \\ ■ \\e is equivalent to the 
standard Euclidean norm on E. 

(3) In terms of the inner product {-,-)e, F\ and E\\ are orthogonal to E—\ and E22, respec- 
tively. 

(4) In terms of the norm \\ ■ \\e, the Lipschitz constant Lp of F with respect to Y satisfies 

L F = 2S!« , (3.9) 
a 

where C2 is as in (HG). 

Note that E 2 is orthogonal to E\ and E = E\ © E 2 . Denote by P and Q (= I — P) the 
projections from E into E\ and E 2 , respectively. Set 



a 

a= 2 



a 5KX\ 



2 a 

Lemma 3.2. (1) For any Y G E 2 , {CY,Y) E < -a\\Y\\ 2 E . 

(2) \\e ct Q\\ E < e~ at for t > 0. 

(3) e ct PY = PY for Y £ E, t > 0. 



(3.10) 
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Proof. (1) and (2) follow from similar arguments as in Lemma 2.3 and Corollary 2.4 in |19] . 
Let us show (3). For Y G E, since PY G E x and f t e ct PY = e ct CPY = 0, we have e ct PY = 
e co pY = py. " □ 

By Lemma 13.21 (2), the constant a in (|3.1Up describes the exponential decay rate of e \qe 
in the new norm. By Lemma 13.11 (4), Lp tends to as a — > oo with respect to the new norm, 
which essentially helps to overcome the difficulty induced from the fact that \i\ — > as a — > 
for i > 1. 

The following lemma will be needed to take care unboundedness of the solutions along the 
direction of the eigenvectors corresponding to pg~. 

Lemma 3.3. Let po = kt]q G E (k is the smallest positive period of g). The random dynamical 
system Y defined in (|3.3p is po-translation invariant in the sense that 

Y{t,uj,Y + po) =Y(t,u,Y ) + po, t>0, wGft, Y €E. 

Proof. Since Cpo = and F(t,co,Y) is po-periodic in Y, Y(t,u,Yo) + po is a solution of (|3.ip 
with initial data Y + Po- Thus, Y(i,w,lo) + Po = Y(t,uj,Y +Po)- □ 

By (|3.3p and Lemma 13.3} (j) is also po-translation invariant. 

Lemma 3.4. For any e > 0, there is tempered random variable r(uj) > such that 

\\z(6 t L))\\ < e e| *!fH for all t G R, uj G fi, (3.11) 

where f{oj) satisfies 

e- eltl f(Lo) < r{Q t u) < e tltl f(io), (el, u G fl (3.12) 

Proof. For j G 1,%, since is a tempered random variable and the mapping t \- > In \ zj(9tui)\ 

is P-a.s. continuous, it follow from Proposition 4.3.3 in [lj that for any ej > there is an tempered 
random variable rj(uj) > such that 

1 

< \zj(uj)\ < rj(u), 



rj(uj) 

where rj(oj) satisfies, for P-a.e. uj € CI, 

e - e ^ rj (uj) < rj(0 t uj) < e^Kjioj), t G R. (3.13) 
Let r(uj) = (rj(uj)) j€li d , u G Cl and take e,- = e, j G 7j d N , then we have 



IKMII < Yl e2e|t|r i(^) = e e| *'lkMII> t € R, wefl. 

Let f{uj) = \\r{uj)\\, uj G fl. Then (l3~TT]l is satisfied and (I3T21) is trivial from (l3~T3|) . □ 



9 



4 Existence of Random Attractor 



In this section, we study the existence of a random attractor. We assume that po = kt]q G E\ 
and 5 G (0, 1] is such that a > 0, where a is as in (|3.10p . We remark in the end of this section 
that such 5 always exists. 

By Lemma [3 .31 and the fact that C has a zero eigenvalue, we will define a random dynamical 
system Y on some cylindrical space induced from the random dynamical system Y on E. Then 
by properties of Y restricted on E2 , we can prove the existence of a global random attractor of 
Y. Thus, we can say that Y has a global random attractor which is unbounded along E\ and 
bounded along E%. Now, we define Y. 

Let T 1 = Ex/pqZ and E = T 1 x E 2 , where p Z = {k Po : k G Z}. For Y G E, let Y := 
Yq (modpo), which is an element of E. Note that, by Lemma l3.3l Y(t, uj, Yb+fcpo) = Y(t,uj,Yo)-\- 
kp , Mk G Z for t > 0, uj G Q and Y G E. With this, we define Y : x S] x E ^ E by setting 

Y(t, w, Y ) = Y(t, u, Y ) (modpo), (4.1) 

where Yo = Yq (modpo). Then Y : R + xl]xE->Eisa random dynamical system. Similarly, 
the random dynamical system <f> defined in (|3.4p also induces a random dynamical system <1? on 
E. By and IfEIjl . * is defined by 

&(t,u},&a)=Y(t,oj,Y ) + z(OtU)) (modpo), t > 0, u G O, (4.2) 

where $0 = 0o {modpo), z(9t<jj) = (0, z(9tUj)) T and Yo = <&o — 2(o;) (modpo). 

Recall that P and Q{= I — P) are the projections from i£ into E\ and i?2; respectively. 

Definition 4.1. Let u> G 12 and 2? : O — > M + be a random variable. A random pseudo-ball 
uj 1 — y B(uj) in E with random radius uj 1— > R(uS) is a set of the form 

u 1 — y B(u) = {b€E: \\Qb\\ E < R(w)}. 

Furthermore, a random set uj h-> B(uj) G E is called pseudo-tempered provided that uj 1— > QB{ui) 
is a tempered random set in E, i.e., for F-a.e.u G fl, 

lim e _/3 * sup{||Q6||_B : b G B{6- t u))} = for all > 0. 

t— >oo 

Clearly, any random pseudo-ball cj 1— > in E has the form w 1— >• E\ X QB(uj), where 

to i-> QB(uj) is a random ball in i?2. Then the measurability of a; 1— > B{u) is trivial. By Definition 
14.11 if uj 1 — y B(u) is a random pseudo-ball in I?, then oj 1— > -B(cj) (modpo) is random bounded 
set in E. And if w h-» -B(cj) is a pseudo-tempered random set in 22, then w h-» (modpo) is 

a tempered random set in E. 

We next show the existence of a global random attractor of the induced random dynamical 
system Y defined in (14. 1|) . 

Theorem 4.2. Xei a > and K > 0. TTien i/ie induced random dynamical system Y defined 
in (|4.ip /10s a global random attractor uj 1— )■ Ao(w). 

Proof. For u; G 12, we obtain from (|3.2p that 

K(t s w,y (w)) = e CT y (^)+ / e^^t^.F^^Fo^)))^. (4.3) 
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The projection of (|4.3p on E2 is 

QY(t, u, Y (u)) = e ct QY (uj) + I e c{t - s) QF(9 s uj, Y(s, uj, Y (uj)))ds. (4.4) 



By replacing uj by 0_tw, it follows from (|4.4p that 

QY(t, 9_ t uj, Y (9_ t oj)) = e ct QY (9_ t uj) + I e c ^QF(9 s ^ t uj, Y(s, 0_ t w, Y {9_ t uj)))ds. 

Jo 

It then follows from Lemma 13.21 and Q 2 = Q that 
\\QY(t,9-tLJ,Y (0-tv))\\E 

< e- at \\QY (6„ t uj)\\E + re- a(t - s) ||F(^_^,y( S ,^,y (^)))kds. ' '" >! 

Jo 

By Lemma 13.41 with e = | and the equivalence of || • \\e and || • || on E, there is a Mi > such 
that 

\\F(e B -tW,Y(s,9-tv,Yo(6-tu)))\\E 

< M 1 ||F(0 a _ t w,Y(s J 0_ t a;,Yo(0-tw 



= Mi(||z(^a;)|| 2 + ||/ - ^(y„(a,e_ t w)) + (1 - a)z(0 s _^)|| 2 ) 5 

< Mx ((3a 2 - 6a + 4)||z(0 s _^)|| 2 + 3||/|| 2 + 3f3 2 c 2 N d ) 1 

< aie^ t_s )f(a;) + a 2 , 



where y u satisfies Y(s, 9-tW, Yo(6~t^)) = (Y u {s,9-tu),Y v (s,9_ t u)) T , a\ = M\ \J 3a 2 — 6a + 4 
and a 2 = M 1 y/3\\f\\ 2 + 3f3 2 c{N d . We then find from (03]) that 



\\QY(t,9- t u,Y (9„ t u;))\\ E < e- at \\QY (9_ t Lo)\\ E + ^(1 - e -f*)fM + ^(1 - e" £ "' ' 

a a 

Now for a; G f2, define 



flo(u>) = r(w) H . 

a a 

Then, for any pseudo-tempered random set ui h-> in £7 and any Yo(^-t w ) 6 B(9_ t u), there 

is a Tb(uj) > such that for i > Tb(uj), 

||QY(t,0_ t w,Y o (0_ t u;))||£ < i?oM, w^, 

which implies 

Y(i,0_ t a;,£(0_ t u;)) C B (u) for all i > T B (w), ui G ft, 

where u; 1— > Bq(uj) is the random pseudo-ball centered at origin with random radius uj h-> Ro(uj). 
Note that cj 1— > Ro{uj) is a random variable since w 1— )■ f(w) is a random variable, then the 
measurability of random pseudo-ball uj h- > Bq(uj) is trivial from Definition 14.11 

For w G fi, let B(w) = (modpo) and Bo(w) = Bq(u;) (modpo), we then have 

Y(t,0_ t w,B(0_ t a;)) C B (w) for all t > T b (uj), uj G O, 
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where Tb(u) = Tb(uj) for oj G f2, i.e., w \-t Bo(w) is the random absorbing set of Y. Moreover, 
uj i — ^ Bo(w) is bounded and closed, hence compact in E, it then follows from Theorem 12.31 that 
Y has a global random attractor u) \— > Aq(w), where 



4>0 T>t 



This completes the proof. 



□ 



Corollary 4.3. Let a > and K > 0. T/ien i/ie induced random dynamical system <& defined 
in (|4.2p /ias a global random attractor u h-> A(w), where A (a;) = Ao(w) + (modpo) for all 



Remark 4.4. fij For any a > and ivT > 0, i/iere is a 5 £ (0, 1] suc/i i/iai a 2 > 25-ftTAi which 
implies a > 0, where a is as in (I3.10P and Ai is i/ie smallest positive eigenvalue of A. 
(2) We say that the random dynamical system Y (or (ft) has a global random attractor in the 
sense that the induced random dynamical system Y (or has a global random attractor, 
and we will say that Y (or (ft) has a global random attractor directly in the sequel. We 
denote the global random attractor ofY and (ft by uj i— >• Aq(uj) and uj i— > A{uS) respectively. 
Indeed, uj \— > Aq(uj) and uj i— > A(uj) satisfy 



Hence a global random attractor ofY (or (ft) is unbounded along the one- dimensional space 
Ei and bounded along the one-codimensional space E^. 

(3) Observe the global attractors of many dissipative systems related to (jl.ip is one- dimensional 
(see |3 \T§ CS UM \M IM El \§Q Similarly, we expect that the random attractor 
uj 1 y A{oS) of (ft is one-dimensional for each provided that a is sufficiently large. We 
prove that this is true in the next section. 

(4) By (2), the system (jl.ip is dissipative along E2 (i.e. it possesses a global random attractor 
which is bounded along E2). In section 6, we will show that (II. ip with sufficiently large 
a and K also has a rotation number and hence all the solutions tend to oscillate with the 
same frequency eventually. 

5 One-dimensional Random Attractor 

In this section, we apply the invariant and inertial manifold theory, in particular, the theory 
established in [7] to show that the random attractor of the random dynamical system (ft generated 
by (jl.ip is one-dimensional (more precisely, is a horizontal curve) provided that a and K are 
sufficiently large (see Remark 14.41 (2) for the random attractor). This method has been applied 
by Chow, Shen and Zhou [8] to systems of coupled first order noisy oscillators and by Shen, Zhou 
and Shen [24] to the stochastic damped sine-Gordon equation. The reader is referred to [31 [5] 
for the theory and application of inertial manifold theory for stochastic evolution equations. 

Assume that po = nrjo and a > 4Lp (see (I3.10P for a and (13. 9p for Lp). Note that the 
condition a > 4Lp indicates that the exponential decay rate of e \qe in the norm || • \\p is 



Proof. It follows from (j4.2j) and Theorem 14.21 



□ 



Ao(w) = Aq(uj) (modpo), A(w) = A{oS) (modpo) 
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larger than four times the Lipschitz constant of F in the norm || ■ It will be seen at the end 
of this section that the condition a > 4L p can be satisfied provided that a and K are sufficiently 
large. 

Definition 5.1. Suppose {^j^eri is a family of maps from E\ to E 2 and n G N. A family of 
graphs oj h-> £{oj) = {{p, <J> w (p)) : p G E±} is said to be a random npo-period horizontal curve if 
oj I— > £(oj) is a random set and {$ w } u gsj satisfy the Lipshitz condition 

I|* w (pi)-* W (P2)||B < \\P1-P2WE for all Pl ,p 2 G E u OJ G O 

and the periodic condition 

^{p + np ) = ^(p) for all p G E 1 , oj G ft. 

Note that for any we!!, £(co) is a deterministic npo-period horizontal curve. When n = 1, 
we simply call it a horizontal curve. 

Lemma 5.2. Xei a > 4Lp. Suppose that uj 1— > £(uj) is a random npQ-period horizontal curve in 
E. Then, uj 1— > Y(t,oj,£(uj)) is also a random npo-period horizontal curve in E for all t > 0. 
Moreover, oj h-> Y(t, 6^ t bJ,£(6-tOj)) is a random npo-period horizontal curve for all t > 0. 

The proof of Lemma 5.2 is similar to that of Lemma 4.2 in [23]. We hence omit it here. 
Choose 7 G (0, §) such that 

«l / 1 1 \ 

< 1, (5.1) 

where is the Lipschitz constant of -F (see ()3.9p ). We remark at the end of this section that 
such a 7 exists provided that a and K are sufficiently large. The main result in this section is 
as follows. 

Theorem 5.3. Assume that a > ALp and there is 7 G (0, |) such that (|5.ip holds. Then the 
global random attractor oj i— )■ Ao(w) 0/ £/ie random dynamical system Y is a random horizontal 
curve. 

Proof. Since equation (|3,ip can be viewed as a deterministic system with a random parameter 
oj G O, we write it here as (|3.ip „, for oj £ £1. We first show that for any fixed oj G 0, (|3.1j) „, has 
a one-dimensional attracting invariant manifold W{oj). 
In order to do so, for fixed well, let 

W(oj) = {Y G E\Y(t,oj,Y ) exists for t < and sup ||e 7 *y(t,o;,yo)lb < °°}- 




We prove that W(oj) is a one-dimensional attracting invariant manifold of (|3,1 
First of all, by the definition of W(oj), it is clear that for any t G R, 

F(i,a;,W(a;)) = W(M> 
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that is, {W (uj)}u£fi is invariant. By the variation of constant formula, Yq £ W(uj) if and only if 
there is Y(t) with Y(Q) = Y , sup t < ||e 7 *y(i)|b < oo, 

Y(t) = e ct i+ f e c(t - s ^PF UJ (s,Y(s))ds+ [ e c ^QF UJ {s, Y{s))ds, i < 0, (5.2) 

JO J -oo 

andY(t,w,y ) = Y(t), where F"{t,Y) = F(9 t u,Y) and £ = P^(O) G E^. For # : (-00, 0] — >• E 
such that sup t<0 ||e 7 *i/(/;)|b < 00, define 

(LH)(t)= f e c{t - s) PH{s)ds+ I e c{t - s) QH{s)ds, t < 0. 



Then 



sup||e 7t (L#)(t)|b < (i + ^_) SU p|| e ^(t)|| B < (± + —L-) sup||e^ff(t)|| Bj 
t<o 7 a ~ 7 t<o \ 7 a - 27 



which means that [|L|| < ^ + a _ 1 2 . Thus, Theorem 3.3 in [7J shows that for any £ G £1, equation 
(|5.2p has a unique solution Y w (i,£) satisfying sup i<0 ||e 7 *y w (/j, £)|b < 00. Let 

/i(u;,£) = QY"(0,Z) = / e- £7a QF w (s,y w (a > 0)<i» ) w G O. 

j —00 

Then, 

W(w)={£ + /i(w,0:£€E£l} 

and is a one dimensional invariant manifold of (|3.ip ,.,. Furthermore, for any e G (0,7), by 

Lemma |3.4| we have 

\\h(9^,0\\E<— r»e £t + -, t>0, (5.3) 
a — e a 



(5.4) 



where ai, 02 is the same as in the proof Theorem 14.21 

To show the attracting property of W(u), we prove for each given u G f2 the existence of a 
stable foliation {W s {Yq,uo) : Yo G VF(u;)} of the invariant manifold W{oj) of (|3.ip „,. Consider 
the following integral equation 

Y (t) = e ct V + j e c ^Q^(s,Y(s) + Y u '( S ^ + h(oj,0)) 

-F w ( s ,y^( s ,e + M^,0)))^ 

+ f e c ^p(F"(s, Y(s) + Y"(s, £ + h(oo, £))) 
J 00 ^ 

-F"(s,Y"( S ,Z + h(oo,0)))ds, t>0, 

where £ + h(u, f) G W(w), 77 = QY(0) G £2 and F w (t, £ + fc(w, £)) := Y(t, 00, £ + 0), * > 
is the solution of (|3.ip with initial data £ + h(u),£) for fixed a; G fi. Theorem 3.4 in [7J shows 
that for any £ € E\ and 77 G E2, equation (j5.4D has a unique solution Y w (i,<!;,?7) satisfying 
sup t > He 7 **™ (£,£,77) lb < 00 and for any £ G 771, 772 G £ 2) 

sup e 7 ' \\Y"(t, £, 771) - £, r? 2 )|b < Af 2 || m - 7/ 2 |b, (5.5) 
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where M 2 = - 2c?w * j ^. Let 

o: V 7 a — / 



^,£,77) =£ + Py" (0,£,77) 

-F"(s,Y"(s,(; + h(u;,Z))))ds. 



Then, £ + £)) = {r/ + /i(w,^) + /i(cj,^,7?) : 77 G £2} is a foliation of at £ + 

Observe that 

Y"(t, e, r/) + r(t, ^ + o) - y w (t, e + o) 

= e a (r? + o + e,»7)-e- M", 0) 
+ ^ e c (*- s ) (>( S) y»(s, t, r?) + y-( s , e + 0)) 

-i^( a ,Y w ( s ,£ + /i( w ,0)))d* 

and 



(5.6) 



= e ct ( V + h(u>, + h(u>, £, 77) - £ - 0) 

+ jf e c (*- s ) (>( S , y^( S , 7/ + + h(u>, £, t/))) 



(5.7) 



Comparing (I5.6P with (I5.7p . we find that 

y w (t,C,7 ? ) = y-(t,7/ + / l ( W ,e) + M^C,7/))-y w (t,C + / l (a;,e)), i>0. (5.8) 

In addition, if 77 = 0, then by the uniqueness of solution of (15. 4|) . y w (i, £, 0) = for t > 0, which 
together with (|5.5p and (|5.8p shows that 



supe^||y^(t,77 + + ^,£,77)) - Y"(t,£ + < M 2 [M| £ (5.9) 

for any £ G i?i and 77 G £?2- Therefore, {W s (Yq,uj) : yo G VT(w)} is a stable foliation of the 
invariant manifold W(u) of (13. ip „, and then W(u) is a one-dimensional attracting invariant 
manifold of (13. lj) „,. 

Next we show that ^4o( w ) = W(u) and Ao(w) is a random horizontal curve. Let u *— > B(ui) be 
any pseudo-tempered random set in E. For any t > and yo G B(0_tu>), there is £(#_tu;,yo) G 
i?i such that 

Y G W 8 (0_ta;, ffl-tU, Y ) + £(0_ t a;, Y ))). 
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Let ri(e-tw) = sup Yo£B{e _ tUj) \\QY - h{9^^{9^ t oj,Y Q ))\\ E . Then by and p, 

\\Y(jt, 9-tV, Y ) - Y(t, 6-tV, W-tOJ, Y ) + £(0_ t a>, Y )))\\ E 

< M 2 e' lt r 1 (e^ t u) 
— > as t — > oo, 

which implies that for oj £ f2, 

d H (Y(t, 9-tfjj, B(6-tfjj)), W{oj)) ^ as t -> oo. 

Therefore, 

= for wGO. 

Moreover, for any random horizontal curve w i-> in E contained in some pseudo-tempered 
random set, 

d H (Y(t,9„ t oj,£(9„ t u>)),A (u;)) as t -> oo 

for every w £ O, which means that lim^oo Y(t, 0_tu;, 1{0^%lS)) C j4o(u;). Since Aq(oj) is one- 
dimensional, we have for oj £ fi, 

4>M = lim y(t,^<a;,^_ t w)). 

t— too 

It then follows from Lemma 15.21 that oj t— >■ Ao(oj) is a random horizontal curve. □ 

Corollary 5.4. Assume that a > ALp and there is 7 6 (0, f ) snc/i £/ia£ (|5.ip holds. Then the 
random attractor oj h- > A(cj) 0/ i/je random dynamical system (f) is a random horizontal curve. 

Proof. It follows from Corollary 14.31 Remark 14.41 and Theorem 15.31 □ 

Remark 5.5. the beginning of this section, we assume that a > <ALp. Since a — — — 1 — — SKXl I 



and Lp = , u>e can tafce a, -fT satisfying ^ 

positive ei 
Note that 



> 8c2 ^ , where Ai is £/je smallest 
positive eigenvalue of A. On the other hand, we need some 7 E (0, |) snc/i i/iai (|5.ip holds. 



a _ 8K\! 
2 a 



2 I 2 



. / 1 1 \ / 1 1 

mm - H — 1 = 1 — h 



(V2 + 1) 



7G(o,|) y 7 a — 27/ I7 a — 27 

which implies that there exist a, K satisfying 

a a 6KX 1 2c 2 |/3|(\/2 + l) 2 8c 2 |/3| 

-- > > . (5.10) 

2 2a a a 

Indeed, let c = 2c2|/3|(v2 + l) 2 , £/zen /or any a > v2c and -K" > i/iere is a 5 > satisfying 

c . . r a 2 - c 1 
KAr <J<mm |^AT' 1 } 

suc/i i/iai fj5. 10f> holds. 
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6 Rotation Number 



In this section, we study the phenomenon of frequency locking, i.e., the existence of a rotation 
number of the coupled second order oscillators with white noises (jl.ip . 



Definition 6.1. The coupled second order system with white noises (jl.ip is said to have a rota- 

T 



tion number p 6 M if, for P-a.e. uj £ Q and each (f>o = (uq, u\) t £ E, the limit lim^oo p f- i '-~ 



exists and 



where r/o is the basis of E\ . 



P(/)(t,UJ,(f)o) 

hm = pjyo, 

t— >00 t 



Note that the rotation number is considered here by restricting <j> on Ei, since 4> is dissipative 
on Ei and limits likewise in Definition 16.11 vanish. From ()3.4p . we have 

P<Kt,u,<M _ PY(t, uj^ojoj)) P(0,z(6 t u;)) T fg ^ 

t t + t 

where <po = (uo, ui) T and Yq(uj) = (uq,ui — z(uj)) t . By Lemma 2.1 in [9j, it is easy to prove that 

lim^oo P(0 ' z( f ta;)) = (0,0) T . Thus, it sufficient to prove the existence of the rotation number 
of the random system (|3.ip . 

Let us show a simple lemma which will be used. For any = s^o € E±, i = 1, 2, we define 

Pi < P2 if si < s 2 . 

Then we have 

Lemma 6.2. Suppose that a > 4Lp. Let £ 6e any deterministic npo-periodic horizontal curve 
(£ satisfies the Lipschitz and periodic condition in Definition \5.1\) . For any Y±, Yi 6 i with 
PYi < PY 2 , there holds 

PY(t,u},Y{) < PY(t,cj,Y 2 ) fort>0,ojen. (6.2) 



The proof of this lemma is similar to that of Lemma 6.3 in |24j . We then omit it here. We 
now have the main result in this section. 

Theorem 6.3. Let a > 4Lp. Then the rotation number of (|3.1|) exists. 

Proof. By the random dynamical system Y defined in (|4.1[) . we define the corresponding skew- 
product semiflow ® t : Q x E — > Q x E for t > by setting 

® t (co,Y ) = (e t uj,Y(t,uj,Y )). 

Obviously, (S7 x E, T x B, (&t)t>o) is a measurable dynamical system, where B = B(E) is the 
Borel a- algebra of E. It also can be verified that there is a measure /i on U x E such that 
(0 x E, T x B, ji, (&t)t>o) becomes an ergodic metric dynamical system (see (6j). Note that 



PY(t,u,Y ) _ PYp 1 
t t t 



f PF(e s u,Y(s,u;,Y ))ds. 
Jo 



17 



Since F(8 s uj, Y(s, u, Ya)+kpo) = F(8 s uj, Y(s, u, Yq)), \/k G Z, we can identify F(6 s u;, Y(s, u), Yq)) 
with F(8 s oj,Y(s, u,Yq)) and write 

F(0 s u;, Y(s,w, Y )) = F(0 s w, Y(a,w, Y )). 

Thus, 

PY(t,co,Y ) PY , 1 '■' 



+ j / PF(^w,Y(a,w,Y ))ds 
* J o 



(6.3) 

PY n i /" 

= — ^ + 7 / F(0 s (^,Y o ))d S , 

where F = P o F G L x (0 x E, J x 6, p). Let t -> oo in ([631), ^r 1 = (0,0) T and by 

Ergodic Theorems in [T], there exist a constant p G R such that 

1 /"* 

lim - / F(& a (oj, Y ))ds = prj , 

t^oo t Jq 

which means 

,. PY(t,cj,Y ) 

hm = pr] . 

t— >oo t 

for [i-a.e.(Lu, Yq) G $7 x E. Thus, there is C SI with P($7*) = 1 such that for any oj G O*, there 
is Y *(w) G F such that 

pr(^,y 1oi)) 

hm — ^ = p?7o- 

t— »oo t 

By Lemma 13.31 we have that for any n G N and well*, 

Py(t,o;,Y *(a;)±npo) PY(t, ^ Y *(a;)) ± np 
hm = hm = pr/Q. (6.4) 

t— Kxi i t— >oo f 

Now for any well* and any 1" G P, there is 71q(w) G N such that 

PY » - n (^)p <PY< PY*(uj) + noHR, 

and there is a ?io(u;)po-periodic horizontal curve lo(u) such that Y *(w) ~~ n o(^)Po> Y, Y *(o;) + 
no(uj)po G lo(u;). Then by Lemma 16.21 we have 

PY(t,u>,Y *(Lj)-n (oo) Po ) < PY(t,u,Y) < PY(t,oo,Y*(w) + n Q (uj)p ), 

which together with (j6.4p implies that for any u G £1* and any Y G F, 

PYft^Y) 
hm = pr] Q . 

t— loo t 

Consequently, for any a.e. and any Y G P, 

PY(t,o,,Y) 
hm = p%. 

t— >oo t 

The theorem is thus proved. □ 

Corollary 6.4. Assume that a > ALp- Then the rotation number of the coupled second order 
system with white noises (jl.ip exists. 

Proof. It follows from (16.11) and Theorem 16.31 □ 
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